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(54) Deformation measuring method and apparatus using electronic speckle pattern 
interferometry 



(57) A deformation measuring method using elec- 
tronic speckle pattern interferometry comprises the 
steps of subtracting an average intensity from the inten- 
sity in a time domain at each point of a speckle pattern 
image so as to compute the cosine component of inten- 
sity; subjecting the cosine component to Hilbert trans- 



form in a temporal domain so as to compute the sine 
component of intensity; determining the arctangent of 
the ratio between thus computed sine and cosine com- 
ponents so as to determine an object phase; carrying 
out an unwrapping operation; and outputting three-di- 
mensional deformation distribution data in a displayable 
mode. 
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Description 

BACKGROUND OF THE INVENTION 
5 Field of the Invention 

[0001] The present invention relates to a deformation measuring method and apparatus using electronic speckle 
pattern interferometry; and, more specifically, to a deformation measuring method and apparatus using electronic 
speckle interferometry which can favorably perform phase unwrapping when analyzing temporal deformation of a dy- 
io namic object 

Description of the Prior Art 

[0002] For measuring the surface shape and deformation of an object in a mirror surface condition by interferometry, 
15 demands for simply determining its phase distribution (corresponding to the surface shape) have been becoming very 
strong especially in the optical and electronic fields due to the advancement of technology in recent years. As techniques 
for determining the phase distribution of such an object to be observed in interferometric measurement in particular, 
those mainly using phase shifting and spatial Fourier transform (Fourier transform introducing therein a spatial carrier; 
ditto in the following) have conventionally been known. 
20 [0003] In general, the phase-shifting method determines a phase distribution of an object according to pattern infor- 
mation of respective interference pattern images, whose phases are shifted from each other by a phase angle which 
is an integral fraction of 2rc, between object light and reference light in an interferometer. The intensity signal in a 
predetermined spatial domain obtained from N interference pattern images equally dividing 2n by N is given by the 
following expression (1): 

25 

W*y) - / 0 [1+7COsW(x,y) + 2ji/7/v}] M>,1,...,(AM) (1) 

where <j>(x,y) is the phase to be determined. On the other hand, !<, is the average intensity, and y is the visibility (mod- 
30 ulation) of interference pattern, both of which are unknown quantities. In a simple example where N = 4, <|>(x,y) is 
determined by the following expression (2): 
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<p(x,y)= tan 




35 <PK x >yj= tan -r —\ (2) 



[0004] The Fourier transform method will now be explained. When one of optical paths of an interferometer is tilted 
by a minute angle 6 about the y axis, a spatial intensity distribution given by the following expression (3) is obtained: 



l(*>y) m M 1 + Y cos {<f>(x,y)+ 2**sin 0/A}] 
- / 0 [l + y cos 2xf Q x}] 

where f 0 = sinGA. is the spatial frequency of carrier fringes generated by tilting the optical path. When the above- 
mentioned expression (3) is Fourier-transformed in the x direction, the following expression (4): 

Af>y)^fj(x 9 y)exp(-j2^x)ix 

(4) 

-'o(/.y)+^(/ + /o>y)+A*(/-/o.y) 

is obtained. 

[0005] The right side of expression (4) indicates that the three terms can be separated from one another on the 
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spatial frequency axis if f 0 is sufficiently large. When f 0 is made sufficiently large so as to take out the second term of 
the right side alone, whereas the other terms are cut by a filter, and the second term is shifted to the origin of the spatial 
frequency domain the following expression (5): 

/,= / oY exp[/<M*,y)] (5) 

is obtained. 

[0006] From the ratio (arctangent) between the real part and imaginary part of the right side of expression (5), the 
10 following expression (6) 



-1 |m LVJ 

*(*.y)=tan (6) 

15 

is obtained. 

[0007] In each of the two conventional methods mentioned above, the unknown quantity <|>(x,y) to be determined is 
obtained independently of the other unknown quantities Iq and y. The phase value is determined as a principal value 
between [-it, n. Fig. 8 shows how the phase is obtained according to expression (6) when the object phase linearly 

20 changes. The denominator and numerator in the above-mentioned expression (6) are cosine and sine functions, re- 
spectively, and the arctangent of their ratio yields a saw-tooth phase distribution in which phase hopping occurs at 
intervals of 2rc. The phase-hopped positions are determined, and the phase value of 2n is added to (or subtracted from) 
data on the right side of each of these positions, so as to correct the phase hopping, whereby a phase distribution in 
proportion to the object shape can be obtained. This processing is known as phase unwrapping. 

25 [0008] In the phase-shifting method, the phase of a point in a two-dimensional interference pattern can be analyzed 
from several intensity data obtained by phase shifting of this point in general. Therefore, this point will not be influenced 
by other points spatially different therefrom. In the Fourier transform method, by contrast, all the data on a line consti- 
tuted by spatial carriers must be subjected to Fourier transform operations, whereby the phase of a point cannot be 
determined independently from other points. 

30 [0009] On the other hand, the object to be observed must stand still throughout a period in which several phase- 
shifted interference pattern images are captured in the phase-shifting method. By contrast, the Fourier transform meth- 
od is considered to be suitable for dynamic phenomena, since phase analysis can be carried out from a single inter- 
ference pattern image only if a carrier component is made beforehand. In other words, the phase-shifting method 
requires the object to be temporally constant though spatial restrictions thereon are loose, whereas the Fourier trans- 

55 form method makes it necessary for the object to have a phase distribution which is spatially moderate enough as 
compared with the spatial carrier period or uniform but allows the object to move temporally. 
[0010] However, these methods are not applicable to cases where the object temporally changes while being spatially 
nonuniform, or changes drastically. For example, the process from plastic deformation of a material to destruction 
thereof is nonlinear, so that a temporal and spatial deformation distribution must be measured when determining dis- 

40 tortions thereof. However, the above-mentioned two methods are theoretically hard to apply to such measurement. 
[0011] Dynamic speckle pattern interferometry has been known as an interferometry method effective in a case 
where object has a rough surface, and such a temporal and spatial fluctuation exists. 

[0012] Speckle pattern interferometry is an interferometric method utilizing a freckle-like pattern (speckle pattern) 
occurring in the observation surface of a rough object illuminated by laser light. In typical imaging systems, the speckle 
45 pattern is considered unfavorable as image noise. However, it carries phase information, so that deformation can be 
estimated from changes in the phase information. In addition, the speckle pattern interferometry enables highly accurate 
deformation measurement with reference to the wavelength of light. 

[0013] Fig. 9 shows a speckle interferometer of a dual illumination type. An object 100 to be observed, which is a 
rough surface object, is illuminated with two luminous fluxes 1 02A, 1 02B from a laser source 1 01 arranged substantially 
50 symmetrical to each other within the x-z plane. The light fields scattered by the object 1 00 form an interference speckle 
pattern on the imaging surface of a CCD camera 103. Thereafter, thus obtained interference speckle pattern image is 
analyzed, whereby phase analysis is carried out to obtain the surface shape of the object 1 00. 
[0014] As the phase analysis technique, subtraction-addition method has been known. 

[0015] In the phase analysis of speckle pattern Images in general, respective speckle patterns before and after the 
55 object 100 is deformed are captured, and differences in intensities of corresponding image points therebetween are 
calculated. The differential intensity approaches zero at places with a stronger correlation, i.e., where the phase change 
caused by deformation is 0 or an integral multiple of 2n, whereas a greater value of differential intensity is obtained at 
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places with a weaker correlation, whereby a correlation pattern corresponding to the amount of deformation can be 
obtained if the absolute value of intensity difference between the two images is calculated. 
[0016] In the above-mentioned subtraction -addition method, information of intensity sum l add is utilized in addition 
to information of intensity difference l sub , whereby a phase is determined independently from the visibility (modulation) y. 
[0017] Here, as shown in Fig. 9, shutters 104A, 104B are placed on the respective optical paths of two luminous 
fluxes 102A, 102B, so as to measure intensity distributions I, (x;t), l 2 (x;t) of individual image in a temporal domain when 
the object is illuminated by only one of the luminous fluxes beforehand. 

[0018] In general, an interference, pattern l(x;t) obtained in an optical system (in-plane deformation or out-of-plane 
deformation system) of a speckle interferometer is represented as in the following expression (7): 

l(x;t) = / 0 (x;f)[1 + Y(x;f)cos(e(x;rM(*;r))] (7) 

where l 0 (x;t) is the average intensity of I, (x;t) and l 2 (x^), 6(x;t) is a speckle random phase, y(x;t) is the modulation, and 
,5 (j)(x^) is the object phase. 

[0019] Subsequently, letting (t - t^ be the intensity before object deformation, and l after (t = t 2 ) be the intensity 
after object deformation, they are respectively represented as the following expressions (8) and (9): 

20 before = f 1 > = 'o t 1 + Y m + *1 )l (8) 

Utter - '(*; fe) - / 0 + y m COS(6 + <fe)] (9) 
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[0020] The difference l sub and sum l add of these intensity patterns are calculated, while local averages are determined 
at the same time, as represented by the following expressions (10) and (11): 

lsub = <l W - 'tetoJ) = c<|sln(0+ A*/2)|X|sin(A<>/2)|>«c1sin(A<|>/2)| (10) 
l add - Waiter + We " 2/ 0» = C<|COS(e+ A*/2)|><|COS(A(|>/2)|> - C|C0S(A<t>/2)| (11) 

where < > indicates a local spatial average, whereas c and c 1 are constants. The 4average intensity l 0 can be determined 
by a local temporal average. A<j> = <j> 2 - ^ and indicates the deformation of object between the time t^ and time t^ 
[0021] The object phase can be determined if the above-mentioned two components are subjected to the arithmetic 
operation represented by the following expression (12): 

A + .2 tan' 1 ^ (12) 

'add 

[0022] When the above-mentioned subtraction-addition method is used, however, the absolute value is taken for 
each of sine and cosine components as represented by the expressions (1 0) and (11 ), whereby the phase curve actually 
obtained for an object phase linearly increasing as shown in Fig. 10 has phase folding points as indicated by a solid 
line in Fig. 10 though not generating the above-mentioned phase hopping. For this matter, each folding section may 
be referred as numeral n, so that phases $ are connected (phase-unwrapped) by using the following expression (13) : 

n - 1, 2, 3, 4,... 
whereby the original object phase $ can be restored. 

[0023] In typical interferometric methods, no absolute values are employed for the sine and cosine components used 
in the arithmetic operation (e.g., the above-mentioned expression (6)) corresponding to the above-mentioned expres- 
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sion (12) for determining the object phase, whereby the resulting phase exhibits saw-tooth phase hopping as shown 
in Fig. 11 (or Fig. 8). Since clear discontinuities occur with a phase jump of 2n> a phase connecting operation for 
restoring the object phase can be executed relatively easily in this case. Namely, it would be more useful if disconti- 
nuities are automatically detected (by software), and an offset of an integral multiple of 2n is added to or subtracted 

5 from this section. In the above-mentioned subtraction-addition method, however, peculiar folding of phase occurs as 
shown in Fig. 10. As a consequence, no discontinuities occur, whereby it is not easy for software to process folded 
sections automatically. In general, speckle noise is remarkable in speckle pattern interferometry as shown in Fig. 12. 
In regions where speckle correlation fringes become dense, i.e., where the amount of deformation is large, the phase 
connecting operation in the subtraction-addition method becomes difficult in particular, whereby automatic processing 

w is likely to impossible in many cases. 

[0024] Also, as can be seen from the above-mentioned expressions (10) and (11), it is fundamentally essentia! for 
the subtraction-addition method to carry out a local spatial average, that degrades the spatial resolution. This makes 
it further difficult to detect phase folding points, whereby the phase connecting operation becomes quite difficult. 
[0025] In view of such difficulty in the phase connecting operation of the subtraction-addition method, the inventors 

15 proposed an improvement (Japanese Unexamined Patent Publication No. 2001-311613) over this technique, which 
has raised the accuracy in detecting phase folding points, but has not been able to detect these poi nts automatically yet. 
[0026] Meanwhile, attention has recently been given to temporal Fourier transform introducing a temporal carrier 
unlike spatial Fourier transform introducing a spatial carrier. If the temporal Fourier transform method is applied to the 
above-mentioned speckle pattern interferometry using a signal having a relatively large noise component, however, 

20 the above-mentioned signal will be hard to process by automatically determining an optimal band-pass filter therefor. 

SUMMARY OF THE INVENTION 

[0027] In view of the circumstances mentioned above, it is an object of the present invention to provide a deformation 
25 measuring method and apparatus which can simplify the detection of phase connecting points in a phase distribution 
curve to such an extent that it is automated when carrying out phase unwrapping in electronic speckle pattern inter- 
ferometry as well, thereby being able to analyze dynamic deformation, vibration, distortion, and the like of an object to 
be observed with a high accuracy. 

[0028] The present invention provides a deformation measuring method using electronic speckle pattern interferom- 
30 etry, the method comprising the steps of determining a phase change curve of a dynamic object to be observed phase- 
wrapped in a predetermined phase range by analysis according to a speckle pattern image carrying phase information 
of the object obtained by using speckle pattern interferometry; and then phase-unwrapping the phase change curve; 

wherein, according to a several of speckle pattern images each obtained at a predetermined time, an intensity 
signal in a temporal domain of each image point is determined, a cosine component of the intensity signal is extracted, 
35 thus extracted cosine component is subjected to Hilbert transform in the temporal domain so as to determine a sine 
component of the intensity signal, and determine a phase change for each image point according to a ratio between 
thus determined sine and cosine components so as to determine a phase change curve of the object. 
[0029] In this case, the Hilbert transform is carried out by using the following expression: 



45 where t and f are times, whereas f(t) is a function of time. 

[0030] Preferably, a phase component eat introducing a temporal carrier for monotonously increasing or decreasing 
a phase term of an intensity signal l(x;t) in a temporal domain for each image point is added to or subtracted from the 
phase term. 

[0031] For example, the phase component cot is added or subtracted by causing a piezoelectric device to move a 
50 mirror surface of a luminous flux reflecting mirror provided for one of two illumination luminous fluxes for generating 

speckle pattern interference so as to change an optical path length of this illumination luminous flux. 

[0032] The deformation measuring method using electronic speckle pattern interferometry in accordance with the 

present invention is particularly useful when the dynamic object is a test piece subjected to a tensile test. 

[0033] The present invention provides a deformation measuring apparatus using electronic speckle pattern interfer- 
es ometry for determining a phase change curve of a dynamic object to be observed phase-wrapped in a predetermined 

phase range by analysis according to a speckle pattern image carrying phase information of the object obtained by 

using speckle pattern interferometry; and then phase-unwrapping the phase change curve; 



40 




the apparatus comprising: 
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intensity signal calculating means for determining an intensity signal l(x;t) in a temporal domain for each image 
point according to a plurality of speckle pattern images each obtained for a predetermined time; 
average component eliminating means for subtracting a predetermined average intensity signal l 0 (x;t) from the 
intensity signal l(x;t) outputted from the intensity signal calculating means so as to compute a cosine component 
5 l c (x;t) of the intensity signal; 

[0034] Hilbert transform operating means for subjecting the cosine component of the intensity signal outputted from 
the average component eliminating means to a Hilbert transform operation in the temporal domain so as to compute 
a sine component of the intensity signal; and 
10 object phase determining means for carrying out an operation for determining an arctangent of a ratio between 

the cosine component of the intensity signal calculated in the average component eliminating means and the sine 
component of the intensity signal calculated in the Hilbert transform operating means so as to determine a phase 
change curve of the object. 

[0035] The apparatus may further comprise temporal carrier superposing means for superposing an appropriate 
is temporal carrier onto the intensity signal l(x;t). For example, the temporal carrier superposing means is optical path 
length difference generating means for generating a predetermined optical path length difference between two illumi- 
nation luminous fluxes for generating speckle pattern interference. 

[0036] For example, the optical path length difference generating means is a luminous flux reflecting mirror configured 
such that a mirror surface thereof is movable with a piezoelectric device. 

20 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0037] 

25 Fig. 1 is a schematic view showing an example of out-of-plane displacement speckle interferometry measuring 

optical system for explaining the concept of the method in accordance with the present invention; 
Fig. 2 is a schematic view showing an example of in-plane displacement speckle interferometry measuring optical 
system for explaining the concept of the method in accordance with the present invention; 
Fig. 3 Is a schematic view showing the configuration of the apparatus in accordance with an embodiment of the 

30 present invention; 

Fig. 4 is a view conceptually showing interference pattern images taken into the apparatus in accordance with the 
above-mentioned embodiment of the present invention; 

Fig. 5 is a block diagram showing individual means included in a computer of the apparatus in accordance with 
the above-mentioned embodiment of the present invention; 
35 Figs. 6A and 6B are charts showing results of analysis of a phase distribution caused by in-plane deformation of 

an object to be observed in Example and Comparative Example, respectively; 

Figs. 7A to 7J are views showing deformation distributions observed by a conventional electronic speckle pattern 
interferometry (ESPI) method in the case where a material such as an aluminum alloy Is pulled; 
Fig. 8 is a graph showing a phase wrapping curve for linearly increasing object phase; 
40 Fig. 9 is a schematic view showing a speckle interferometer (measurement optical system) of a dual luminous flux 

irradiation type; 

Fig. 10 is a chart showing a phase folding state in the subtraction-addition method; 
Fig. 11 is a chart showing a phase folding state in a typical interferometry method; and 
Fig. 12 is a view showing speckle noise in an image obtained by speckle interferometry. 

45 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 

[0038] In the following, an embodiment of the deformation measuring method using electronic speckle pattern inter- 
ferometry in accordance with the present invention will be explained with reference to the drawings. 

50 

Basic Configuration 

[0039] Figs. 1 and 2 show respective examples of optical systems for measuring out-of-plane and in-plane displace- 
ments of an object to be observed by using speckle pattern interferometry. 
55 [0040] In the interference optical system apparatus for measuring out-of-plane displacement shown in Fig. 1 , laser 
light emitted from a semiconductor laser light source (LD) 11 is partly scattered by a diffuser 13 arranged in front of an 
object 1 2 to be observed, so as to form a speckle pattern. On the other hand, the part of laser light transmitted through 
the diffuser 13 illuminates the object 12 and is scattered by the surface thereof, so as to form a speckle pattern in a 
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manner similar to that mentioned above. These two scattered light components form an image of a random interference 
speckle pattern on the imaging surface of a CCD camera 1 4, which photoelectrical^ reads out this interference pattern. 
In the interference optical system apparatus for measuring in-plane displacement shown in Fig. 2, on the other hand, 
the object 12 is illuminated with two laser light beams substantially symmetrica! to each other about a normal of the 
5 object which are split by a half mirror 17, so that speckle patterns independent from each other are formed, whereby 
an image of an interference speckle pattern is formed on the imaging surface of a CCD camera 1 4 due to the interference 
of these two speckle patterns. Thereafter, this interference speckle pattern is photoelectrically read out in a manner 
similar to that mentioned above. PZT 15 will be explained later. 

[0041] In each of the apparatus mentioned above, the intensity l(x;t) of the random interference pattern detected by 
10 the CCD camera 14 is represented as the following expression (14): 

/(x;f) = / 0 (x;r)[1 + Y (x;r)cos(e(x;f) + 4>(x;f))] (14) 

15 where x is a coordinate, t is the time, l 0 (x;t) is the average intensity, -y(x;t) is the degree of modulation, 6(x;t) is the 
random speckle phase, and <|>(x;t) is the object phase accompanying the displacement of the object to be observed. 
[0042] This embodiment aims at determining the object phase <|>(x;t) with a high accuracy. 

[0043] Therefore, in this embodiment, the average intensity l 0 (x;t) is initially determined by local temporal averaging 
with respect to the intensity in the above-mentioned expression (14), and then is subtracted from the interference 
20 pattern intensity l(x;t), so as to take out only the cosine component of interference pattern intensity. Letting l c (x;t) be 
thus taken-out component, it can be represented as the following expression (15): 
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/ c (x;r) = /(x;f) - / 0 (x) = /(x;rK/(x;r)>, <* cos(e(x;f)+(|>(x;f)) (15) 

[0044] Subsequently, the temporal signal represented by the above-mentioned expression (1 5) is subjected to Hilbert 
transform. 

[0045] The Hilbert transform with respect to a function f (t) is defined by the following expression (1 6): 

F(r) = i/X.[/(0] = -iL^' (16) 

[0046] By subjecting the cosine component of the above-mentioned expression (15) to Hilbert transform, the sine 
component represented by the following expression (17): 

/ s (x;f) = HT,[/ C (x;f)] « sin (G(x;f) + 4>(x;f)) (17) 

is obtained. 

[0047] The ratio between the above-mentioned cosine component l c (x;t) and sine component yxjt), and the arctan- 
gent of this ratio is determined, whereby the object phase is determined as the following expression (18): 



0(x;t) + <l>(x;t)<= tan" 1 



(18) 



[0048] Normally, the change in speckle random phase 6(x;t) is slower than the change in object phase, whereby it 
can be approximated as in the following expression (19): 

6(x;f) s 6 (x;0) (19) 
[0049] Therefore, the target object phase is determined as in the following expression (20): 
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tan" 1 

5 

[0050] When the observation time ^ is so long that expression (19) does not hold, the total observation period ^ is 
divided into n sections with a shorter period ^ (to = ny in which the speckle phase can be considered constant, and 
the object phase is determined in each section by substituting 6(x;0) with efxiky, k = 0, 1 , .... n-1 . Namely, the reference 
phase is updated as appropriate. Finally, the individual sections are connected together, whereby a dynamic phenom- 
10 enon over a long period of time can be analyzed. 

[0051] Thus determined phase curve in the temporal domain is in a state whose phase value is wrapped (folded) 
within -n and % as principal values as mentioned above. Therefore, thus wrapped (folded) phase curve is subjected to 
phase-unwrapping for restoring it to the original curve. 

[0052] By carrying out the foregoing arithmetic operation for each image point on the speckle pattern interference 
is image, deformation of the object to be observed can be analyzed by quantifying it dynamically with a high accuracy. 

Determination of Polarity of Phase Term 

[0053] Meanwhile, in the explanation of principle mentioned above, uncertainty occurs in the polarity of phase term 
20 as can be seen from expression (15). For eliminating this uncertainty, a known phase component art is intentionally 
introduced into the phase term of expression (15). Namely, a phase component art determined such that 6(x;t)+<l>(x;t) 
+ert, which is the phase term in the following expression (21), monotonously increases or decreases is introduced as 
so-called temporal carrier component 

/(^;0= / O (x;f)[1+7(x;0cos(8(x;0+<|)(x;f)+(of)] (21) 

[0054] Specifically, it can be carried out by various known techniques, for example, such as one modulating the 
current injected to the semiconductor laser light source 11 as shown in Fig. 1 and one changing the optical path length 
30 by moving the mirror 1 6 with the piezoelectric transducer (PZT) 1 5 as shown in Fig. 2. 

Apparatus Configuration 

[0055] Rg. 3 is a schematic diagram showing a deformation measuring apparatus using electronic speckle pattern 
35 interferometry in accordance with an embodiment of the present invention. This apparatus is an example of apparatus 
for carrying out in-plane displacement measurement of an object to be observed. 

[0056] As explained in Rg. 2, the object 12 to be observed is illuminated with two laser light beams, substantially 
symmetrical to each other with respect to a normal of the object, split by the half mirror 1 7, whereby random speckle 
patterns are formed independently from each other due to scattering of the two laser light beams from the object 12. 

40 The interference between these two speckle patterns forms an interference pattern on the imaging surface of the CCD 
camera 14. Here, the piezoelectric transducer (PZT) 15 is used for changing the optical path length of one of the 
systems of the interferometer in order to introduce the above-mentioned temporal carrier component The piezoelectric 
transducer 1 5, CCD camera 1 4, and semiconductor laser light source (LD) 1 1 are under the control of a timing controller 
21 , whereby the image data of interference pattern captured by the CCD camera 1 4 is stored into a memory within a 

45 computer 22. In the computer 22, the image data is subjected to phase analysis processing using Hilbert transform, 
so as to obtain information Items such as the amount of distortion of the object. Thus obtained information items con- 
cerning the object are shown on a display 23. 

[0057] The speckle pattern images are captured at predetermined intervals as shown in Fig. 4, whereby the temporal 
change of intensity l(x,y;t) is obtained at each image point within the images (see the graph at the right end of Rg. 4). 
so [0058] Here, the intensity l(x,y;t) at one point (x,y) within the images is represented by the following expression (22): 

/(*,y;0 = tt^yiOtl+Yt^y^Jcoste^yiO+^Jx.yiO+fflf)] (22) 

55 where art is the temporal carrier component; © =» 0 when this component is not introduced. 

[0059] As shown in the block diagram of Rg. 5, the arithmetic operations in the computer 22 are sequentially carried 
out in an average component eliminating means 32, a Hilbert transform operating means 33, an object phase deriving 



'ten 



-e(x,o) 



(20) 
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means 34, and a phase connecting means 35 with respect to the intensity l(x,y;t) (hereinafter simply referred to as I 
(x;t)) as the input signal. 

[0060] Here, the intensity l(x;t) is one in which an appropriate carrier component is introduced thereon by temporal 
carrier introducing means such as: PZT mentioned above. 
5 [0061] The inputted intensity l(x;t) is initially subjected to an arithmetic operation for eliminating the average intensity 
l 0 (x;t) by using the above-mentioned expression (15) in the average component eliminating means 32. 
[0062] The Hilbert transform operating means 33 computes the sine component by Hilbert transform using the above- 
mentioned expression (17). 

[0063] The object phase deriving means 34 carries out an arithmetic operation for determining the arctangent of the 
10 ratio between the cosine and sine components by using the above-mentioned expression (18), and then determines 
the object phase by using the above-mentioned expression (20). 

[0064] Thereafter, the phase connecting means 35 carries out an unwrapping operation for the object phase curve, 
thereby determining dynamic deformation information of the object. Namely, as mentioned above, phase-hopping po- 
sitions of the phase curve are determined, and a phase value of 2rc is added to (or subtracted from) the data on the 
15 right side of each of thus determined positions, so as to correct the phase hopping, thereby yielding a phase distribution 
in proportion to the deformation of object. 

[0065] Finally, a three-dimensional deformation distribution forming means 36 yields a three-dimensional deformation 
distribution such as the one shown in Fig. 6A from the phase curve obtained after the phase unwrapping operation for 
each image point, and its data is outputted after being converted into a mode which can be represented on the display 

20 23 three-dimensionally. 

[0066] The above-mentioned means 32 to 36 are constructed by programs stored in the computer 22. 
[0067] Since the method of the present invention is particularly effective when employed in deformation measurement 
for carrying out a tensile test (or fatigue test) of a material, for example, effects obtained when the method of the present 
invention is applied to a process of analyzing temporally changing deformation of an aluminum alloy in a tensile test 

25 thereof will now be explained with reference to specific examples. 

Examples 

[0068] In general, when a material such as an aluminum alloy is pulled, uniform stretches occur in the elastic defor- 
30 mation area thereof. When they are observed by electronic speckle pattern interferometry (ESPI), uniform fringes such 
as those shown in Figs. 7A to 7C are seen. The direction and density of fringes vary depending on the position due to 
the fact that rotations are caused by unexpected installation errors of the test piece and the like in addition to the 
deformation. If the deformation exceeds the yield point, thereby reaching the plastic deformation area, the deformation 
will be localized in a narrow region tilted by about 45° with respect to the pulling direction, which is known as slip band, 
35 in general. When it is observed by electronic speckle pattern interferometry (ESPI), fringes are concentrated in the 
localized part as shown in Figs. 7D to 7F. 

[0069] Further, when observed with a relatively long differential time for the rate of deformation, a band-like part in 
which, an inner fringe structure seems to have disappeared emerges and moves up and down at a substantially constant 
speed. When a sufficiently short differential time is taken for this part, a fringe structure appears. When this part is 
40 observed in detail under magnification, it can be seen that local deformation develops nonlinearly in a complicated 
manner as shown in Figs. 7G to 7J (observed at intervals of 3.6 seconds from Fig. 7G to 7J). Such a conventional 
speckle pattern interferometry method (ESPI) can qualitatively grasp temporal changes of deformation but cannot 
quantitatively determine the magnitude of deformation. 

[0070] When the method of the present invention is applied to the processing for analyzing the deformation, by 
45 contrast, the magnitude of deformation can quantitatively be grasped with a high accuracy in a simple manner even 
when the S/N ratio is relatively low. 

[0071] Using the apparatus of the above-mentioned embodiment, a phase distribution caused by in-plane deforma- 
tion of an object to be observed was analyzed under the following experimental condition. Fig. 6A shows the results. 
[0072] The results obtained by using the above-mentioned subtraction-addition method under the same condition 
50 are shown in Fig. 6B as a comparative example. 

(1) Experimental Condition 

Object: aluminum alloy (with a thickness of 5 mm) 
55 Pulling rate: 0.5 nm/sec 

Observation area: 80 mm x 30 mm 

Frame rate: 30 frames/sec 

Illumination light: semiconductor laser light 
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Illumination light wavelength: 532 nm 
(2) Experimental Results 

As can be seen when Figs. 6A and 6B are compared with each other, the subtraction-addition method yielded 
5 quite erroneous analysis data since the phase connection could not be carried out correctly, whereas the magnitude 

of deformation could be obtained accurately when the method of the present invention was used. 

Comparison with Cases Using Other Techniques 

w [0073] As mentioned above, various problems occur when other techniques are used in place of the method of the 
present invention in the case where deformation of an object to be observed in a temporal domain is measured by 
using electronic speckle pattern interferometry. 

[0074] Namely, it is difficult for phase-shifting or spatial Fourier transform method to handle a speckle pattern which 
changes temporally and spatially. On the other hand, it is difficult for temporal Fourier transform method to automatically 
15 determine an optimal band-pass filter for a signal having a low S/N ratio in speckle pattern interferometry. It is difficult 
for the subtraction-addition method to find out the phase folding points, so that automatic processing is impossible in 
many cases, whereby only quite erroneous analysis results can be obtained as shown in the above-mentioned com- 
parative example except for cases with favorable data acquiring conditions. 

[0075] When the method of the present invention is used, by contrast, favorable analysis results can be obtained 
20 without the fear of generating problems such as those of the above-mentioned conventional techniques. 

Modification of Embodiment 

[0076] The method in accordance with the embodiment of the present invention can be modified in various manners. 

25 For example, unwrapped phase data may be subjected to smoothing by use of a low-pass filter or median filtering. 
[0077] The present invention can be applied to dynamic objects in general. For example, it is applicable not only to 
fatigue tests of materials in general including items such as the above-mentioned tensile test of materials, but also to 
minute dynamic changes of various constituent parts of animals and plants, such as human tissues in particular. 
[0078] As explained in detail in the foregoing, the deformation measuring method and apparatus using electronic 

30 speckle pattern interferometry in accordance with the present invention extracts the cosine component of the interfer- 
ence speckle pattern intensity of the object to be observed, subjects thus extracted cosine component to Hilbert trans- 
form in a time domain so as to determine the sine component of the intensity, and determine the phase at each image 
point according to thus determined sine and cosine components, thereby yielding a phase distribution curve of the 
object. Therefore, a phase distribution curve can easily be determined even in an object which changes temporally 

35 and spatially, while phase connecting points can easily be specified even from a signal including a large amount of 
noise typical in speckle pattern interferometry, so that automation is possible, whereby quite highly accurate analysis 
results can be obtained. 

[0079] A deformation measuring method using electronic speckle pattern interferometry comprises the steps of sub- 
tracting an average intensity from the intensity in a time domain at each point of a speckle pattern image so as to 
40 compute the cosine component of intensity; subjecting the cosine component to Hilbert transform in a temporal domain 
so as to compute the sine component of intensity; determining the arctangent of the ratio between thus computed sine 
and cosine components so as to determine an object phase; carrying out an unwrapping operation; and outputting 
three-dimensional deformation distribution data in a displayable mode. 

45 

Claims 

1. A deformation measuring method using electronic speckle pattern interferometry, said method comprising the 
steps of determining a phase change curve of a dynamic object to be observed phase-wrapped in a predetermined 

50 phase range by analysis according to a speckle pattern image carrying phase information of said object obtained 

by using electronic speckle pattern interferometry; and then phase-unwrapping said phase change curve; 

wherein, according to a plurality of speckle pattern images each obtained at a predetermined time, an intensity 
signal in a temporal domain of each image point is determined, a cosine component of said intensity signal is 
extracted, thus extracted cosine component is subjected to Hilbert transform in said temporal domain so as to • 

55 determine a sine component of said intensity signal, and determine a phase change for each image point according 

to a ratio between thus determined sine and cosine components so as to determine a phase change curve of said 
object. 
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A deformation measuring method according to claim 1 , wherein said Hibert deformation is carried out by using the 
following expression: 



where t and V are times, whereas f(t) is a function of time. 

A deformation measuring method according to claim 1 or 2, wherein a phase component tot introducing a temporal 
carrier for monotonously increasing or decreasing a phase term of an intensity signal l(x;t) in a temporal domain 
for each image point is added to or subtracted from said phase term. 

A deformation measuring method according to claim 3, wherein said phase component <ot is added or subtracted 
by causing a piezoelectric transducer to move a mirror surface of a luminous flux reflecting mirror provided for one 
of two illumination luminous fluxes for generating speckle pattern interference so as to change an optical path 
length of said one illumination luminous flux. 

A deformation measuring method according to one of claims 1 to 4, 
wherein said dynamic object is a test piece subjected to a tensile test. 

A deformation measuring apparatus using electronic speckle pattern interferometry for determining a phase change 
curve of a dynamic object to be observed phase-wrapped in a predetermined phase range by analysis according 
to a speckle pattern image carrying phase information of said object obtained by using speckle pattern interfer- 
ometry; and then phase-unwrapping said phase change curve; 
said apparatus comprising: 

intensity signal calculating means for determining an intensity signal l(x;t) in a temporal domain for each image 
point according to a plurality of speckle pattern images each obtained for a predetermined time; 
average component eliminating means for subtracting a predetermined average intensity signal l 0 (x;t) from 
said intensity signal l(x;t) outputted from said intensity signal calculating means so as to compute a cosine 
component l c (x;t) of said intensity signal; 

Hilbert transform operating means for subjecting said cosine component of said intensity signal outputted from 
said average component eliminating means to a Hilbert transform operation in said temporal domain so as to 
compute a sine component of said intensity signal; and 

object phase determining means for carrying out an operation for determining an arctangent of a ratio between 
said cosine component of said intensity signal calculated in said average component eliminating means and 
said sine component of said intensity signal calculated in said Hilbert transform operating means so as to 
determine a phase change curve of said object. 

A deformation measuring apparatus according to claim 6, further comprising temporal carrier superposing means 
for superposing an appropriate temporal carrier onto the intensity signal l(x;t), said temporal carrier superposing 
means comprising optical path length difference generating means for generating a predetermined optical path 
length difference between two illumination luminous fluxes for generating speckle pattern interference. 

A deformation measuring apparatus according to claim 7, wherein said optical path length difference generating 
means is a luminous flux reflecting mirror configured such that a mirror surface thereof is movable with a piezoe- 
lectric device. 
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